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I. INTRODUCTION 


The National Transonic Facility (NTF) is a complex wind tunnel facility (Figure 1 A.). The 
purpose of this user’s guide is to briefly describe the facility, the data system, the instrumentation 
used to acquire research data and to define some of the acronyms commonly used at the facility. 
The computational methods and equations will be discussed in detail and references will be listed 
for the user who needs additional technical information. The customer should find this guide 
helpful both in preparing for a test and during the operational phase of the test. It is intended that 
this document answer many of the general questions about how research data is acquired and pro- 
cessed in the National Transonic Facility. 

II. FACILITY DESCRIPTION 

The NTF is a fan-driven, closed-circuit, continuous-flow pressurized wind tunnel (Figure IB). 
The test section is approximately 8.2 X 8.2 feet square and 25 feet long with longitudinal slots in 
the top and bottom walls (6 slots each, 6 percent open). Test section configuration variables are 
top and bottom wall divergence, reentry flap angle, and initial diffuser angle (i.e. model support 
wall angle). Refer to Figures 2A-2L for tunnel operating envelopes. 

The NTF test gas may be air or nitrogen. For the air mode of operation, cooling is obtained by 
water-fed cooling coils located at the upstream end of the settling chamber. The cooling capacity 
of this system is 3 1 MW. For the nitrogen mode of operation, heat is removed by evaporating liq- 
uid nitrogen which is sprayed into the tunnel circuit upstream of the fan. By using liquid nitrogen 
as a coolant, the tunnel test temperature range is variable from -250 to 150 °F (1 16 to 340 °K ). 
The maximum cooling capacity at -250 ° F is 115 MW. Thermal insulation has been installed 
internal to the pressure shell to minimize energy consumption. 

To facilitate rapid model access for configuration modifications during a test series with a min- 
imum loss of nitrogen and time, model access tubes are provided (Figure 3). The tubes are 
inserted from either side of the plenum and close around the model to form a thermally insulated 
enclosure. This enclosure can be conditioned for safe personnel entry while the plenum and tun- 
nel circuit remain cold. 

The tunnel drive system consists of two variable speed induction motors with a combined 
power output of 35 MW (50 MW short term overload), a two speed gear box, and a synchronous 
motor with a power output of 3 1 MW (45 MW short term overload). The compressor consists of 
a fixed pitch, single stage fan with variable inlet guide vanes. The drive system can be operated in 
either of two modes: (1) the synchronous mode where the fan is operated at the constant synchro- 
nous speed of 360 RPM and the inlet guide vanes are varied to achieve the desired compression 
ratio ; and (2) the variable fan speed mode where the induction motor(s) are used to drive the fan 
over a range of rotational speeds (up to 600 RPM) to achieve the desired compression ratio. Dif- 
ferent combinations of fan speed, inlet guide vane angle and gear ratio are used by the tunnel 
operator depending on tunnel pressure and temperature. 
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III. DATA SYSTEM DESCRIPTION 


The NTF operation is centered around energy efficiency, high productivity and safety. It is 
unique in the aeronautics community because of its ability to produce high quality data at flight 
Reynolds Numbers. The NTF data system plays an important role in achieving the high produc- 
tivity and high quality data which is so important to our customers, both NASA and industry. 
Timely dissemination of data to our customers is a goal of this facility, and a reliable and efficient 
data system is necessary to achieve this goal. 

The National Transonic Facility Open Architecture Data System (Figures 4A-4C) consists of 
three data acquisition/reduction systems in the Real-Time Systems (RTS) group: (1) the Research 
Computer System (RCS), the Process Computer System (PCS) and the Model Preparation Area 
(MPA); (2) a Host Computer System (HCS); and (3) sixteen Networked Display System Con- 
soles (NDSC). Together, these computer systems provide the functionality for acquiring, process- 
ing, recording and displaying on-line and off-line data from the model, the tunnel and the model 
preparation areas. 

The RTS group is based on the MODCOMP REAL/STAR 2000 computer running a real-time 
operating system called REAL/IX. This operating system is based on the standard Unix (AT&T 
System V) with real-time extensions and provides a user friendly environment which is modular, 
scalable and has distributed system components. The use of this open computer architecture 
allows many commercially available, “off the shelf’ hardware and software products to be used 
for system expansions and upgrades. The CPU functions are handled by four 25-MHz Motorola 
88100 CPU’s in a tightly coupled, fully symmetrical processor configuration. 

The PCS data acquisition subsystem is responsible for collecting the data from several differ- 
ent data acquisition units (DAU). The primary analog inputs consist of a Neff 620/500 controller 
and a Neff 620/600 analog input subsystem with up to 256 analog channels and 16 digital chan- 
nels. Other analog inputs come over a serial RS-232 interface from a PC-based Temperature Mon- 
itor Scanner system which monitors up to 640 temperatures for tunnel safety and structural 
integrity assessment. Three sequencers, consisting of Programmable Logic Controllers (PLC) 
which control and sequence tunnel interlocks, transmit up to 4096 discrete digital state inputs to 
the PCS. The PCS also receives data over a high-speed SCRAMNet reflective memory system. 
This data comes from three control microprocessors and from the RCS system and is used to con- 
trol model attitude. 

Since safety is a major concern in a high pressure cryogenic environment, the PCS has an 
elaborate alarm system. This system provides a channel-based alarm processing capability with 
characteristics which are configurable in the test setup file to provide an independent configura- 
tion for each channel. When a channel enters or exits an alarm state, a message is logged to the 
alarm display window, the system log window and the system log files. A channel can also be 
configured to activate visual and/or audible annunciators. These annunciators are located on the 
Test Engineer’s panel and their on/off status is controlled by the alarm system software. Each 
channel has a deadband parameter which reduces unwanted alarms, without affecting safety. 
Optionally, a response procedure can be displayed for corrective action. 
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The PCS system is divided into eight different functional areas as shown in the PCS block dia- 
gram (Figure 5A). The functional areas are initialization, test setup, system control, calibration, 
data collection, computation, alarming and data display. Each is a separate and distinct area of 
PCS but necessary for the operational success of the total system. 

The RCS (Figure 5B) data acquisition subsystem is responsible for collecting the data from a 
similar Neff DAU containing 128 analog channels and 16 digital channels. This system continu- 
ously acquires, computes, performs alarm checks and displays model data and tunnel parameters. 
The RCS computer also receives model and tunnel wall pressure data from an electronically 
scanned pressure system (ESP 8400) across an IEEE 488 interface. Standard force and moment 
data are measured with strain gage balances which are calibrated to account for first and second 
order interactions. Other measurements may be made such as pressures from the model, tunnel 
walls, boundary-layer rakes, internal flow rakes, and wing gages. Model angles are measured 
with either an on-board accelerometer package or the accelerometer located in the arcsector. 
Model roll angle is measured by either a resolver or potentiometer. Since NTF is a cryogenic 
wind tunnel, temperatures are monitored in many regions of the tunnel and on the model. The 
RCS computer also processes data from thermocouples, platinum resistance thermometers ( PRT) 
and various other temperature measuring devices. The subsystems used for computing Mach 
number will be described in a later section (III A). During operation, two tasks run in the RCS 
computer: (1) TCYCL acquires, averages, reduces and displays cyclic or “real-time” data and (2) 
TPOINT acquires, averages, reduces, records and displays static (point-based) data. The RCS 
system has the same eight functional areas as the PCS system and the block diagram is the same 
with the exception of the input sources. 

To minimize the time required for model installation and instrumentation verification, one of 
the three model preparation areas can be used for pretest setup. The MPA data acquisition sub- 
system is patterned after the RCS and runs a subset of the RCS software. The computer system is 
used for accelerometer calibrations, balance loadings, sting bending calculations and pressure ori- 
fice leak checks. Following model assembly and ambient temperature checkout, cryogenic mod- 
els undergo further checkout in the cryogenic chamber located in MPA 3. The data acquired 
during the cryogenic checkout is later used for temperature corrections during the actual tunnel 
test. The computer system used in the NTF model bays is very effective in model pretest prepara- 
tion, instrumentation verification and software checkout. 

As stated above, each of the RTS computers interfaces with the Neff 620/500 controller and the 
Neff 620/600 DAU to acquire analog and digital data. The Neff 620/600 is a high-speed unit that 
measures multi-point, low-level and high-level analog signals. This unit provides the system with 
the following features and functions (Reference 14): 

Analog to digital conversion 
16-bit resolution (including sign) 

Accuracy of 0.02% of full scale ±2 microvolts at constant temperature after 
autocalibration 

- Programmable gain post amps with steps 1,2,4,8,16 and 32 

- Autoranging capability with throughput of 50 kHz (currently used on fixed range) 
Maximum throughput of 100 kHz 
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Automatic zero and full-scale calibration 
Precision internal voltage calibration source 

Alphanumeric display for viewing channel values on-line and off-line. 

The Neff 600 unit uses the Neff 620650 amplifier/filter circuit cards to signal condition the 
incoming data signals before they are digitized. Each card contains four channels and has the fol- 
lowing features and functions: 

Programmable gain preamps with steps 1,8,64, and 512 

Programmable filter selection, using 4-pole Butterworth low-pass filters, with steps 
of 1,10,100 and 1000 Hz 

Full scale voltage ranges from 5 millivolts to 10 volts 
Input protection of 100 VDC or peak AC 

- Common range rejection ratio of 66 dB to 120 dB depending on gain and filter settings 
Buffered wide band output available for each channel for external use. 

As stated above, one of the capabilities of the Neff 620/600 is an automatic internal self-cali- 
bration for all data channels (defined in the acquisition setup file) using all gain and filter steps 
available. Zero and full scale corrections are determined for each channel. This information is 
stored in the Neff internal calibration memory, written to calibration disk files and written to test 
recording file if enabled (raw data file). Both the number of scans and the averaging period is 
configurable. At the present time, the usual scan rate for a discrete data point is 10 scans/sec and 
the averaging period is one second. 

Although the purpose of this document is to provide information about the research data sys- 
tem, information has also been included about the PCS because tunnel related information is 
often used by both the test engineer and the research engineer. Since these computers are on a 
local area network (LAN), information from the PCS and the RCS can be passed back and forth 
and also mixed on displays and plots. This is helpful particularly during cryogenic and/or high 
pressure operations since the tunnel environment affects the quality of the research data. The out- 
puts from pressure gages and thermocouples located throughout the tunnel are recorded on the 
PCS system and can be displayed along with the model data from RCS if desired. 


IV. TUNNEL PARAMETERS 


A. Parameter Measurements 

The NTF research Mach system was designed to measure tunnel total and static pressures to a 
high degree of accuracy such that the resulting Mach Number accuracy would be ±0.002 over the 
range from 0.2 to 1.2. The transducers for this system are fused quartz absolute gages (RUSKA) 
and are staged to achieve the required accuracy levels. For static pressure measurement there are 
five gages with the following full scale values: 150, 100, 50, 30, 15 psia. For total pressure there 
are only four gages with the 15 psia gage omitted since the tunnel is not designed to operate at 
total pressures less than one atmosphere (14.7 psia). Solenoid valves provide automatic over- 
pressurization protection at 90% of full scale. The data system monitors the position of these 
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valves to pick the pressure values from the most sensitive gages. 

The gages for this system have a laboratory quoted accuracy of±0.012 % of reading plus 
±0.006 % of full scale. The sensitivities of the transducers are determined by off-site laboratory 
calibrations. Two point in-place calibrations are made daily with these calibrations consisting of a 
near vacuum point ( < 200 microns) to determine and adjust the gage offset value and a baromet- 
ric pressure point that is used in combination with the vacuum point as a check on sensitivity 
deviation. 

A correction is made to the total pressure measurement as a result of the transducers and the 
probe orifice not being located on the centerline of the test section. The orifice is located near the 
bottom of the tunnel in the settling chamber and the gages are at a lower elevation. This correc- 
tion is of the form: 

P T,CL = ^ T,Meas~ (^7^2 + ^l^l) 

where 

Xj = specific weight of gas at tunnel test temperature and pressure, f — 

V/ 

\ 2 - specific weight of gas at ambient temperature and tunnel pressure, 
h x = 20.58/if , height of column of gas at tunnel temperature and pressure 

h 2 = 2.0ft , height of column of gas at ambient temperature and tunnel pressure. 

See sketch below for details. 
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An additional pressure measurement system (Mach system) is currently being verified and will 
be used as a redundant system. The flow reference system (FRS) is based on PS I 8400 pressure 
system technology. This system consists of four pressure gages (two absolute and two differen- 
tial). The absolute gages are resonant quartz transducers (Digiquartz) and have a full scale range 
of 150 PSI. These transducers are accurate to ±0.015 % of full scale. The two differential pres- 
sure gages are piezo-resistor micro-machined silicon sensors (HASS) and have full scale ranges 
of +30 and +15 PSID respectively. The accuracy of these transducers is also ±0.015 % of full 
scale. One of the absolute gages measures tunnel total pressure, ( P T ) from a pressure probe 
mounted on the wall of the contraction cone at the horizontal centerline of the tunnel. This instal- 
lation removes the necessity for an elevation correction to the P T measured by this system. The 

other absolute gage measures the tunnel reference static pressure, ( P s ) . Both of the delta pres- 
sure gages measure P T -P S . The most sensitive A P gage (within range) is used in combination 
with P T to calculate Mach Number. If neither of the A P gages is within range, then the pressures 

from the two absolute gages are used in the calculation. It is anticipated that this system will be as 
accurate as the Ruska based system now in use. The gages for this system are housed in modules 
that are sent to the laboratory to be calibrated on a semi-annual basis. This system may eventually 
be used as the primary Mach system and may also be used for tunnel Mach control. 

The prime measurement of tunnel stagnation temperature is made with PRT probe mounted in 
the the reservoir section of the tunnel near the screens. Currently, there are several of these 1000 
ohm probes, two of which are used at any given time. These probes have an accuracy of approx- 
imately ±0.05 “/f and a response time of 1 to 2 seconds. 

The data reduction program utilizes the Callender- Van Dusen equation to compute temperature 
(Reference 1): 

3- xa-itt-.x*)') 

where 

a, (3, 8 are coefficients in the laboratory calibration 

R 0 is the resistance at 0°C 

R t is the resistance at temperature, T . 


The tunnel PRT instrumentation system consists of a constant current power supply for the 
PRT and a high accuracy calibration resistor ( R c = 1000 ohms ). The PRT and calibration resistor 
are connected in series as shown on the next page. 
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With this electrical arrangement, V T and V c are measured and -^= f 7/) x because 

An V r ) V An / 


V T V 0 V c 

1 D _ V „ J D _ C 


^ 7 ^= — , /? 0 = — , and R c = — . With the resistance ratio determined, the stagnation temperature 

is calculated from the Callender- Van Dusen equation which is iteratively solved in the data reduc- 
tion software. 


Measurements of tunnel stagnation temperature are also made with a copper-constantan ther- 
mocouple which has an accuracy of ±1°AT at cryogenic temperatures. The thermocouple has a 
response time of 1 to 2 seconds. Copper-constantan thermocouples are also used to measure 
model, balance, and tunnel component temperatures. All these thermocouples are referenced to 
the ice point (0°C) and the millivolt outputs are converted to temperature by using polynomial 
fits of the form 


T= A 0 + (AjV) +(a 2 V 2 ] + [a 3 V 3 ) + (a 4 V 4 J 

where V is in millivolts and T is in °K . Each polynomial covers a specific millivolt range as indi- 
cated in the table below. 
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COEFFICIENTS FOR COPPER-CONSTANTAN THERMOCOUPLES 
{Reference Temperature, 273. 15°^ or 32 °F ) 


MV 

RANGE 

-6.059 < -5.587 

-5.588 <-3.574 

-3.575 < 2.864 

2.865 < 20.684 

AO 

4020.6970000 

97.65544600 

273.15000000 

274.00707000 

A1 

1441.7082000 

-143.83491000 

25.85383400 

25.18987600 

A2 

43.3593140 

-62.52461700 

-0.71469546 

-0.54212165 

A3 

-34.5629920 

-10.01594900 

-0.07198166 

0.01734062 

A4 

-3.3589826 

-0.63851732 

-0.01351040 

-0.00025029 


B. Calculation of Hinnel Conditions 

Gas Properties 

In the NTF data reduction software, the gaseous properties for both air and nitrogen are math- 
ematically represented by the Beattie-Bridgeman thermal equation of state which has the form 


= r «m-en mi 

V 2 J IV 2 


(Reference 2) 


and by ideal-gas (i.e. zero pressure) specific heats of the form 


C p G, G 2 G^ 9 _3 Gop 2 ^ 

r r (/-l) 


where p= — (Reference 3). 


The constants for these two descriptions are given in the following tables. The first table contains 
the equations of state (Beattie -Bridgeman) constants. 


*= ' 4 o( 1 
8 = B 0 ( 1 
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where 



NITROGEN 

AIR 

1 

1.34450 

1.30120 

a 

0.02617 

0.01931 

B 0 

0.05046 

0.04611 

m 

-0.00691 

-0.01101 

c 

4.2 X 10 4 

4.34 X 10 4 

w 

28.01340 

28.96000 


R (gas constant) = 0.0820539 liter-atm / mole - K 
V(specific volume)= 1/p where p is density 
W = molecular weight. 

The above constants apply when the quantities have the following units: P(atm), T(K), and 
V(liters/mole). 

The second table contains the ideal gas specific heat constants. 



NITROGEN 

AIR 

Gi 

-7.352104012 X 10 2 

0.0 

g 2 

3.422399804 X 10 1 

0.0 

G 3 

-5.576482846 X 10' 1 

0.0 

g 4 

3.504042283 

3.5 

G 5 

-1.733901851 X 10' 5 

0.0 

G 6 

1.746508498 X 10' 8 

0.0 

g 7 

-3.568920335 X 10' 12 

0.0 

OO 

O 

1.005387228 

0.0 

G 9 

3.353406100 X 10 3 

0.0 
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In addition to the equation of state and specific heat representations of the two gases, the NTF 
software also mathmatically characterizes the viscosity of the gases. 

For nitrogen, the equation is p= p 0 (T) + p e (p,7) (Reference 4) 

where 
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M7)= X 2 , 7 ’ 0 " 3 * 

/ = 1 

7 

IV 

is 1 


Zi 

7.41653 X 10 1 

x, 

2.30835 X 10' 1 

Z 2 

-1.58344 

X 2 

-9.36362 X 10' 1 

Z 3 

3.85308 X 10' 3 

X 3 

9.033392 

Z 4 

8.01337 X 10‘ 4 

X 4 

-4.18321 X 10 1 

Z 5 

-8.92031 X 10' 7 

x 5 

1.08976 X 10 2 

Z 6 

8.90597 X 10' 10 

x 6 

-1.29139 X 10 1 

Z 7 

-5.37794 X 10' 13 

X 7 

5.97821 X 10 1 

Z 8 

1.73983 X 10' 16 



z 9 

-2.30840 X 10' 20 




These constants apply when the quantities have the following units :p (10 3 gm/cm-s), T (K),and 
p (gm/cm 3 ). 


For air, Sutherland’s equation is 


F = 


1.45820xl0" 6 


j3/2 

(7+110.33) 



Thermodynamic Equations 

The thermodynamic property equations for entropy, enthalpy, specific heats and speed of sound 
are as follows: 
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S{p,T) = 



dT- [/? x In {pRT) ] 



dp 


rT C 


,o\ 


where ^ r 0 ,Amb + | r [ / j ^ is an arbitrary constant = 25.5308 (jou!es/mole)and is represented 
by G,, in the data reduction software. 


« (P. * T) = ♦ f o [( A) - 0 )( !)J ♦ X C > 


r T 0 

where H T Amb + C dT is an arbitrary constant = 14.9124 (joules/mole) and is represented by 

p J T 0 

G ]0 in the data reduction software. 


Cv(p, T) = 




dpi 


37 2 ) 


\dp 


Cp(p,T) = 


C(p,T) = 





J 


where 


Amb subscript 
0 superscript 
T 0 subscript 


implies p = 1 atmosphere 
implies p = 0 

implies a constant temperature of 298.15 °K . 


Since air is only used at near ambient temperature in the NTF, a constant ideal gas specific 
heat representation is adequate. Previous research (Reference 5) has shown that for the operating 
temperature and pressure range of NTF, this reasonably simple representation of the properties of 
nitrogen is in close agreement with much more complex representations that were designed for 
high accuracy over a much broader pressure and temperature range. In addition, prior research 
(Reference 6 and 7) has shown that various types of flows (isentropic, normal shock, and bound- 
ary layer) are not significantly different from those same flows calculated for an ideal diatomic 
gas. While the nondimensional flow parameters that characterize these flows are essentially the 
same, proving that cryogenic nitrogen is a good test medium, it is necessary to retain this more 
complex representation of the gas properties to calculate the various dimensional flow quantities 
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such as density, speed of sound, and dynamic pressure. 
Procedure 


With measurements of free-stream total pressure, total temperature and a static reference pres- 
sure, the other important reference conditions such as Mach Number, Reynolds Number, and 
dynamic pressure can be calculated by using basic thermodynamic equations (Reference 8) which 
express properties such as entropy, enthalpy, specific heat, and speed of sound in terms of the 
equation of state, its derivatives and the ideal gas specific heat values. These equations along with 
the assumption that the flow is isentropic (constant entropy) and adiabatic (no heat transfer) make 

possible an iterative solution of the adiabatic energy equation H T = H + ( V*/2) where H is 

enthalpy and V is velocity. Once this solution is complete, all of the desired reference quantities 
can be calculated. The various steps in the calculation procedure are as follows: 

1 . Input the measured test parameters and the gas type (air or nitrogen). 

a. P { , stagnation pressure 

b. P s , plenum static pressure 

c. T t , stagnation temperature 

2. Calculate other stagnation properties. 

a. p = p l (P l , T t ) .density 

b. S= S t (p,, T t ) , entropy 

c. H= H t ( P t , T t , S t ) , enthalpy 

3. Iteratively solve for static temperature assuming isentropic flow. 

a. Approximate static temperature, T, using ideal-gas expansion from P T to P, T. 

b. Calculate static density, p= p (P, T) 

c. Calculate static entropy, S- S (p, T) 

d. Obtain incremental variation in entropy with temperature, AS/ AT, by incrementing T 

e. Check for convergence, S= S t ± ^ S t x ](f 8 J 

(1) If convergence criterion has not been met, adjust T by T= T+ ( (5 - S) / ( AS/ AT) ) 

and return to step b. 

(2) If convergence criterion has been met, static temperature, T, and density, p , are 
now known. 

4. Determine other tunnel parameters. 

a. H- H(P,p,T) , enthalpy 

b. C— C (p, T) , speed of sound 

c. V= [2 {H t -H) l/2 J , flow velocity 

d. M- V/C, Mach number 
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e. 4= p (p, T) , viscosity 

f. R e - (pW)/p, Reynolds number 

0 

g. q= (pV)/2, dynamic pressure 

C. Calibration Corrections 

A correction is made to the reference Mach Number (Mrpp) based on the results of centerline 
static pipe calibration tests in the NTF. The correction is a function of M REF and Reynolds 
Number per foot ( REYNFT). A double interpolation of the following table values is performed 
and the correction, A M , is added to to obtain the free-stream Mach Number ( ). 


AM CALIBRATION CORRECTION 


m ref 

REYNFT 

6 X 10 6 

14 X 10 6 

30 X 10 6 

60 X 10 6 

100 X 10 6 

BOX 10 6 

0.0 

0.0 X 10' 3 

0.0 X 10' 3 

0.0 X 10' 3 

0.0 X 10' 3 

0.0 X 10' 3 

0.0 X 10' 3 


-1.3 

-1.3 

-1.3 

-1.3 

-1.3 

-1.3 

0.30 

-1.7 

-1.9 

-2.8 

-2.8 

-3.5 

-4.1 

0.40 

-2.1 

-2.1 

-2.2 

-2.2 

-2.3 

-2.3 

0.50 

-2.3 

-2.4 

-2.6 


-3.5 

-3.9 


-2.7 

-2.7 

-2.8 

-3.0 

-3.3 

-3.5 

0.70 

-2.5 

-2.6 

-2.8 

-3.1 

-3.6 

-3.9 

0.80 

-2.5 

-2.7 

-3.0 

-3.6 : 

-4.4 

-5.0 

0.90 

-2.8 

-3.0 

-3.4 

-4.1 

-5.0 

-5.7 

0.95 

-2.5 

-4.3 

-7.6 

-4.7 

-4.7 

-4.7 

0.98 

-1.9 

-4.5 

-7.7 

-7.8 

- 7.8 

- 7.8 


-2.3 

-5.1 

-6.0 

- 7.8 

-7.3 

-6.7 

1.02 

-2.6 

-3.7 


- 7.5 

- 7.6 


1.05 

-2.5 

-5.8 

-6.6 

-8.3 

-8.3 

-8.3 


-2.1 

-4.7 

-5.7 

-7.5 

- 7.5 

-7.5 

1.15 

-8.4 

-8.4 

-8.4 

-8.4 

-8.4 

-8.4 
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Once AM has been determined, can be calculated as follows: 


M ref 


+ AM 


Another correction which is based on the NTF static pipe calibration data is the buoyancy cor- 
rection ( A CD buoy ) to the drag coefficient (CD ). This correction is also a function of Mref and 

REYNFT. A double interpolation of the following table values is performed to obtain — . For 

dX 

this calculation it is assumed that the top and bottom test section wall angles are zero, the re-entry 
flap angles are zero, and the initial diffuser angles are -1.76°. 


dM/dX CALIBRATION CORRECTIONS 


m ref 

REYNFT 

6X 10 6 

14 X 10 6 

30 X 10 6 

60 X 10 6 

100 X 10 6 

BOX 10 6 

0.0 

0.0 X 10‘ 4 

0.0 X 10' 4 

0.0 X 10' 4 

0.0 x 10- 4 

0.0 X 10' 4 

0.0 X 10' 4 

0.20 

-0.989 

-0.885 

-0.677 

-0.287 

0.233 

0.623 

0.30 

-1.436 

-1.371 

-1.241 

-0.996 

-0.670 

-0.426 

0.40 

-1.849 

-1.771 

-1.613 

-1.317 

-0.923 

-0.628 

0.50 

-1.992 

-1.914 

-1.757 

-1.462 

-1.069 

-0.775 

0.60 

-1.851 

-1.778 

-1.634 

-1.363 

-1.001 

-0.730 

0.70 

-1.506 

-1.446 

-1.326 

-1.100 

-0.800 

-0.575 


-0.675 

-0.628 

-0.534 

-0.358 

-0.123 

0.053 

0.90 

0.184 

0.150 

0.080 

-0.049 

^0.222 

-0.352 

0.95 

0.636 

-0.550 

-2.050 

-0.800 

-0.800 

-0.800 

0.98 

1.250 

-0.300 

-4.000 

-4.400 

-4.400 

-4.400 

1.00 

1.500 

-0.200 

-3.000 

-5.100 

-4.200 

-3.350 

1.02 

-0.900 

-0.350 

-0.900 

-2.200 

-2.200 

-2.200 

1.05 

3.900 

-0.280 

-0.500 

-0.900 

-0.900 

-0.900 

1.10 

12.400 

7.800 

7.600 

7.200 

7.200 

7.200 

1.15 

7.940 

7.940 

7.940 

7.940 

7.940 

7.940 
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Using ^ from the table, the buoyancy correction and the corrected drag coefficient can then 
be computed using the following equations: 



where 


CD t — CD , + A CD l. 

corrected uncorrected buoy 

Vol = model volume, ft 3 
S ref = model reference area, ft 2 
= free-stream Mach Number 

^22 = tunnel longitudinal Mach gradient, 1/ft. 


This correction is applied directly to the drag coefficient instead of the axial force so that it does 
not affect any of the other coefficients through Euler rotations. 


V. MODEL PARAMETERS 
A. Model Orientation 

Critical measurements in any wind-tunnel model test are those of model attitude with respect to 
the gravity axis system. These measurements, along with a knowledge of the tunnel flow direc- 
tion with respect to gravity, makes possible a determination of the model orientation with respect 
to the wind axis. Existing techniques at Langley Research Center involve the use of accelerome- 
ters to sense the attitude of the model (Reference 9). For the cryogenic environment of NTF, it is 
necessary to provide heated containers for the accelerometers so that they remain at a constant 
temperature ( « 160 °F). There are many considerations when choosing the proper on-board 
package to be used at NTF. The space and wiring requirements of a package must be considered. 

Another inertial device used in NTF is the electrolytic bubble which detects the angle between 
its base and local horizontal. This device when installed in an on-board accelerometer package is 
used to establish “wind-off zero attitude” during the test. 

The NTF acceleromenter packages NTF10-X and NTF12-X measure elevation angle of the 
model. In addition, the NTF12-X package contains two electrolytic bubbles measuring zero ele- 
vation both upright and inverted. Since the NTF10-X has no bubbles, it should only be used when 
the space available prohibits the the use of the “12” package. The output of the accelerometer is a 
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function of the sine of the angle so it has its best accuracy at small angles. The accelerometers 
will measure within ±0.01° for small angles. 

Prior to a wind-tunnel test, the model support system is assembled in one of the three model 
preparation areas. Sting and balance deflection constants are determined for normal, side, pitch 
and yaw loadings. This also provides a good check on the balance and the data system. If the 
model is to be cryogenically tested, it will either be assembled in or moved to MPA 3 where the 
model/sting combination will be thermally cycled to cold test temperatures. Once all calibrations 
are completed in the MPA, the model/sting assembly can be transported to the test section by a 
model handling cart. The balance is again checked in the tunnel by hand loading. 

Once a model has been installed in the tunnel, a reference plate is attached to the model which 
establishes the model reference axes. The AOA package should then be powered and heated to 
160°F for 30 minutes before proceeding with the calibration. The pitch and roll 0° output 
(model level) for the AOA device is recorded. If an NTF-12 package is being used, the zeroing 
pot on the bubble readout is adjusted for the upright bubble channel. A second calibration point is 
recorded at pitch =0° and roll= 90° . The leveling plate should have provisions for inverted lev- 
eling of the model so that a third point can be recorded at pitch= 0° and roll= 1 80° . At this 
point, the zeroing pot on the bubble will be adjusted for the inverted bubble channel. Now all the 
necessary information has been obtained for calculating the calibration constants. The following 
equations are used for this calculation: 


where 

Eleva 

E 

Bias 

Sens 

Omax 

Az 

Roll 


Eleva- asin j ~ Omax x sin (Az + Roll) 

= elevation angle (usually referred to as “pitch”angle in software) 
= output of accelerometer ( volts) 

= electrical bias of accelerometer (volts) 

= sensitivity of accelerometer (volts/g) 

= misalignment of accelerometer (deg) 

= roll orientation of misalignment (deg) 

= roll attitude of model (deg) 


The sensitivity of the accelerometer is acquired from the laboratory calibration. Bias, Omax 
and Az are computed from the in-place tunnel measurements described above. The equations are 
as follows: 


Bia (E ° + £ '» 0> 

2 

where 

^lso = voltage output at 180° roll. 
E q = voltage output at 0° roll and 

Az and Omax are then calculated using this bias. 
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Az= atan 


and 


where 


" asin ( (E 0 - Bias ) / Sens ) 
asin ( (Eqq - Bias) /Sens) 


Omax= 


asin ( (E 0 - Bias) / Sens) 
sin ( Az ) 


E 90 = voltage output at 90° roll. 


Basic three dimensional model support is provided by an aft-mounted sting. This sting is 
attached to a vertically mounted arcsector which is driven by a hydraulic cylinder. The roll mech- 
anism provides the interface between arcsector and sting (Reference 10). The pitch range of the 
the arcsector is from - 11 ° to 19° , but sting knuckles can be used to obtain higher pitch angles. 
The roll range is -180° to 180° , and sideslip angles are achieved by using pitch/roll combina- 
tions. Computer software is available to compute the pitch (0) and roll (<b) angles given input 
angles of a (model angle of attack) and (5 , (model angle of sideslip). Design load capabilities of 
the model support system are discussed in Reference 10. 

The measured angles with respect to the gravity axis that are obtained from the primary orien- 
tation system are converted to angles with respect to the other axis systems by applying a series of 
Euler rotations (APPENDIX A and APPENDIX B) The model accelerometer package is always 
considered the primary system but a secondary system consisting of an accelerometer in the arc- 
sector can be selected by using the angle definition code. See Figures 6 A through 6 E for descrip- 
tions of how these two systems can be used in the software. 

B. Model Forces and Moments 

In order to measure three dimensional model forces and moments in the NTF, a series of special 
balances ( APPENDIX C) have been designed with metals that have acceptable strength and frac- 
ture toughness at cryogenic temperatures. These moment-type balances have two sets of strain 
gages mounted near both ends of a beam section. The sum of the gage outputs is proportional to 
the force (normal and side) and the difference in the outputs is proportional to the moment (pitch 
and yaw). There are individual gage outputs for axial force and rolling moment. The gages for 
these balances were carefully selected, matched and temperature compensated such that sensitiv- 
ity and no-load output variations are in the order of 1 % of full scale output for the temperature 
range 300 to 100° K. The changes in sensitivity with temperature are measured in a laboratory 
cryogenic chamber and sensitivity shift constants are calculated for the balance. The zero offset 
changes with temperature are measured here at the NTF in the cryogenic chamber of MPA 3. The 
chamber is cycled from ambient temperature to approximately 100 °K . The sensitivity and zero- 
load voltage is recorded (usually at five minute intervals) during the MPA cryogenic cycle. The 
data are fit with a second order curve and coefficients are calculated for each gage. These coeffi- 
cients are included in the standard balance data reduction to further reduce errors due to tempera- 
ture variations. For ambient temperature testing, the data reduction software is equipped to 
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handle both force and moment type balances. 

The temperature correction philosophy indicates that all balance readings (wind-off zero, 
weight tare data, and wind-on data) that are taken at a given temperature should be corrected to a 
reading that would have been obtained at a fixed ambient temperature. This fixed temperature is 
the temperature at which the balance is calibrated (295 ° K). The data for each gage is correlated 
with balance temperature as measured by the temperature measuring device nearest the gage. All 
NTF balances have three temperature devices located at the front, middle and rear of the balance. 
The actual placement of these devices is shown in Figure 7. The zero shift due to temperature 
level is characterized by the following equation 


* v ez = <VC,7, + C 2 J? 

where: 

C 0 , Cj, C 2 = temperature compensation coefficients 

7, = temperature from temperature device for i th gage (1-6) 

EZ = subscript denoting electrical zero. 

The effect of sensitivity shift on load voltage is expressed as 

V L = ^(Z.,295) X 

where: 

5C,= 1+55,(7,-295) 

55, = sensitivity change per degree change in temperature 
= (5-5 295 )/(5 295 xA7) 

L = subscript denoting load. 

Thus, the equation for the measured bridge voltage, V T is: 

V T = V (EZ, 295) + A V EZ + ( V Tare + ^ Aero ) L 

where: 


V (ez, 295 ) = bridge electrical zero at T = 295 ° K 
v t.„ = voltage due to model weight 
V Aero ~ v °ltage due to aerodynamic loads. 

Since V Tare and V Aero in the above equation are voltages due to load, they are related to their 
calibration temperature values as follows: 
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V Tare + V Aero= ( V 7V,r, + V AeJ (Z .,295) X 5C 

The temperature compensation algorithm philosophy is to correct the measured voltage, V T , to 
what would have been measured at the calibration temperature (295° K). 


V(T, 295)- V (EZ, 295) + ^(EZ, 295) + ( V T.re +V A era ) L /(SC) 


where 

AV^z 295 ) approaches zero at the calibration temperature. 

F~ ( ( V / , 7 .. Q ,,) , . , , ) x Sens, ( ( V. TOO c,) , , „)xSens 

i v v (T.295)' (wind-on)' ' vv (T,295y (wind-off) ’ ' 

where Sens . = balance primary sensitivities. 

Even though the data reduction software applies temperature compensation to each of the bal- 
ance components, test procedures call for taking wind-off zeroes at test conditions. This mini- 
mizes the effect of the zero offset correction. 

The result of the above procedure is a set of balance forces and moments that have been cor- 
rected to a temperature reference of 295 ° K. The forces and moments must now be corrected for 
1 st and 2 nd order interactions, for model weight tares, and other corrections such as base and 
chamber pressure. Also, these forces and moments must be transferred to the model and wind 
axis systems and then non-dimensionalized in coefficient form. Although the details of these pro- 
cedures are very similar for all wind tunnels, the balance data reduction equations including all 
corrections are given in APPENDIX A. 


C. Model Pressures 

Steady-state pressure measurements in NTF are made with several different types of instru- 
mentation depending on the accuracy requrements. For multi-channel applications such as mea- 
surements on the test section walls and on models with multiple pressure orifices, ESP 
instrumentation is used. 

The ESP system consists of several key components (Figure 8 ). The multi-transducer modules 
(available with 16, 32, 48, and 64 ports) are housed in thermal enclosures in the cryogenic envi- 
ronment of NTF and the temperature is maintained at near ambient levels. The modules are inter- 
faced with the system processor which contains the scanner digitizing units (SDU) through the 
combination of the remote power supply (S81FC) and the scanner junction unit (SJU). The small 
size of the SJU allows the transducer voltage to be controlled at or very near the model. This unit 
is also housed in a thermal enclosure. The modules are pneumatically connected to the model ori- 
fices, tunnel reference pressure, the local processor with its slot mounted pressure calibrate units 
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(PCU) and the bias control unit. This unit biases the control line pressure (shifts module to run/ 
calibrate position) to valves that are approximately 100 PSI above the tunnel reference pressure. 

At NTF, this system is configured such that the host computer commands the ESP system to 
read/send raw voltage at a prescribed scan rate over a specified period of time. This period of 
time (nominally, 1 second) is the same acquisition period used for other types of data such as bal- 
ance data and tunnel parameter data. The maximum scan rate is dependant on the number of 
ports on each module in the system. 

The average voltage from each module is converted to engineering units (PSID) using the fol- 
lowing 4 th order polynomial equation 


A P= C 0 + C, V+ C 2 V 2 + cy + C 4 V 4 


where 


A P= 


P - P 

orifice ref 


V = average voltage from a given port/module 
Cj = calibration coefficients 

P re r- module reference pressure 


The calibration coefficients are determined from two types of calibrations. Full calibrations 
involve the PCU sending out five pressures (0, -FS, -.5FS, +.5FS, FS) where FS indicates the full 
scale pressure of a module. At certain tunnel conditions, the negative pressures as indicated can- 
not be achieved. In these cases, there is an algorithm used to predict what the substitute pressure 
should be . The purpose of this algorithm is to calculate the value of five reasonably spaced cali- 
bration pressures when the tunnel pressure drops below the critical pressure for that size module. 
Then the calibration pressures are calculated internally in the software and operator intervention 
is not required. 

The algorithm is as follows: 

(1) Calculate the critical pressure for each size module: 

P crit = ( P ref~ P miJ wher e P min = minimum pressure of the vacuum system (P min ~ 5 psi) 

(2) Calculate pressure interval for module calibration based on comparison of P to FS 
pressure of a module: 
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If 

(Pcri,>FS) 

then 

P 

min 

= P ref-FS 

and 

P 

max 

= P ref 

If 

(Pcm<FS) 

then 

p 

min 

= P 

mtn,v 

and 

p 

max 

= P ref 


(3) Set calibration pressures based on P cnl : 


If 

( P cn,* FS ) 

then: 

PI 

= 0 


P2 

= -FS 


P3 

= -FS / 2 


P4 

= +FS / 2 


P5 

= +FS 


If 


then: 

PI 

= 0 


P2 

= (p . -p 

v min 

ref) 

P3 

= (P -p 

v min 

ref) ^ 

P4 

= + FS / 2 


P5 

= +FS 



FS. 

FS 


These calibrations take a minimum of 25 seconds per pressure (2.1 minutes total time) and are 
usually only done at wind-off conditions. Zero calibrations involve the PCU sending a zero A P 
across the transducer. The higher order coefficients are retained from the last full calibration and 
the offset coefficient (C 0 ) is updated to satisfy the above equation. During calibration , on-line 
diagnostics are used to indicate that the system(s) is performing correctly and within expected tol- 
erances. The various diagnostics are given in APPENDIX D. In order to check the ESP sys- 
tem^) during tunnel operation, at least one port on each module is pneumatically connected such 
that A P will equal zero for that port. When the indicated pressure on that reference port exceeds 
the accuracy tolerance (±0.19 % FS), a zero calibration is usually performed. 

Tunnel wall pressures (16 rows, 352 “good” orifices) can be measured and used to make wall 
interference assessments. The orifice layout and the standard module hook-up remain essentially 
the same from test to test. 

Pressure data may be presented in coefficient form. The coefficients are calculated by the 
following equation: 
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where 


4 , 


and 


P„ = free-stream static pressure 
= free-stream dynamic pressure. 


VI. OPTIONAL DATA CORRECTIONS AND CALCULATIONS 
A. Sting-bending Temperature Corrections 

An investigation was conducted at the NTF to document the change in sting-balance deflec- 
tions from ambient to cryogenic temperatures (Reference 15). A method of reducing the required 
calibration time is to obtain only ambient temperature sting- balance bending data and to correct 
for changes in the modulus of elasticity with temperature. The following equation is used to cor- 
rect the sting bending constants: 


K. = K t (A 0 + A.T r ) 

* mYn v 1 a 


where: 


Kj= yaw deflection due to side force 
K- 2 = yaw deflection due to yawing moment 
K 3 = pitch deflection due to normal force 
K 4 = pitch deflection due to pitching moment 
K 5 = roll deflection due to rolling moment 
K^= yaw deflection due to rolling moment 
K 7 = roll deflection due to side force 
Kg= roll deflection due to yawing moment 

CRYO indicates cryogenic temperatures 

AMB indicates ambient temperature 

T r = rear balance thermocouple temperature in °F. 

Using the above equation over the temperature range from 75 °F to -250 °F, a decrease of 
approximately 3.7% in bending can be expected. For the most prevalent sting material, Vasco- 
max 200, A 0 = 0.990746 and A, = 0.0001322. If other materials are used, these constants can 
be adjusted in the data reduction set-up file. 



B. Ttonnel Flow Angularity Correction 

Model orientation can be corrected for tunnel flow angularity in the pitch (upflow) direction 
and/or the yaw (crossflow) direction. Refer to Appendix B for a discussion of the Euler rotation 
equations. Upflow angle is determined by measuring normal force coefficient, C N with the 

model upright and inverted at the same tunnel conditions. The inverted polar angle of attack, a 
values are multiplied by (-1) and a first order curve fit is performed on both polars over a specified 
angle range. This angle range is chosen based on the model. The intersection point of these two 
curves becomes the upflow angle, The method is shown below. 


^~'N upr Qupr U/>* UPR 
C N, nv = C O INV +C \ wv a iNV 

To determine flow angularity, C N 


(1st order fit) 
(1st order fit) 
= C N and a 

/v INV 


a INV , therefore: 


and 


+ C i ctf — C n + C , 


a. 


a f = 




( C i, 



where 

C Q = offset 
C, = slope 

UPR indicates upright 
INV indicates inverted 
CLj = upflow angle. 


Currently, the Data Analysis System (DAS) software (Reference 13) is used to compute flow 
angularity. The data is first processed using the off-line data reduction software. A file is created 
containing the values to be used for computing the flow angle correction. Once the upflow angle 
has been determined another pass is made throught the data reduction software to correct for tun- 
nel flow angularity. The present method is somewhat awkward and time consuming and a one 
pass data reduction method is planned for the near future. This method will be self-contained and 
will not require the use of the DAS software. 

In the NTF, there is no mechanism for conducting an angle sweep in the yaw direction, there- 
fore, only two points (<|> = 90° and <p = -90° ) can be obtained to determine crossflow. This nor- 
mally is not done since it would not give an accurate representation of tunnel crossflow angle. 
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C. Mass Flow Calculation 


Mass flow ratio calculations for flow-through nacelles are often made at the NTF. First the 
average pressure coefficient is calculated from n static pressure orifices in the nacelle. 


CPIAV = 


(CP l + CP 2 +...CP n ) 

n 


Using CPIAV, the average internal Mach Number, MIAV can be calculated. 




MIAV = 


CPIAV xq +P 

l oo < 


- 1 


00 / 


where 

P T = free-stream total pressure 
q „ = free-stream dynamic pressure 

P^ = free-stream static pressure 


The mass flow ratio (nacelle flow/free-stream flow), MFR is then calculated using the follow- 
ing equation: 


MFR = KVx 


MIA V 

J [ MIAV* y 


M 


r M 2 '* 
\ 00 

1+ T 


) 



where 

M^ = free-stream Mach Number 

ATF = nacelle flow area 

AHL = free-stream reference area 

KV = viscous correction parameter (non-dimensional) based on chord Reynolds Number. 


24 



D. Calculation of Model Deformation Parameter 


The amount that a model deforms is predicted by the ratio of free-stream dynamic pressure, q M 
to the modulus of elasticity, E, of the model metal. Aeroelastic studies in the NTF can be per- 
formed by holding the Mach Number and Reynolds Number constant and varying the q m /E 

parameter. Studies can also be made with a constant model shape (constant q^/E) while varying 

Reynolds Number. The modulus of elasticity (Young’s modulus) for a particular model material 
is a function of temperature, therefore, this parameter is calculated by the following equation: 

E = 144.0 [E 0 + (£, x T t ) + [ E 2 x T* t J + ^£ 3 x J] 

= free-stream dynamic pressure in PSF 
T t = tunnel total temperature in °F. 

For the two model materials most frequently used in the NTF, the coefficients are as follows: 

(1) Vascomax 

E 0 = 26.47 X 10 6 
E, = -3468.0 
E 2 = 0.0 

e 3 = o.o 

(2) A286 

E 0 = 28.371 X 10 6 
Ej = 754.9 
E 2 =22.218 
E 3 = .025944 


VII. DATA MANAGEMENT 
A. Data Acquisition and Reduction 

The RCS must be configured for each test program based on the requirements documented in 
Administrative Instruction Procedure-29 (AIP-29) and the Test Plan developed by the Test Engi- 
neer. Once the configuration file has been generated, the system is started through the menu and 
the process of defining the required configuration environment will automatically be accom- 
plished. One set-up file contains all the necessary information for the acquisition, recording , and 
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on-line data reduction. This file defines the types of instrumentation to be used, the DAU chan- 
nels which have been selected, the ESP configuration and in some cases, calibration information 
for various tunnel devices. It also contains retrieval information for balance calibrations and tem- 
perature compensation constants. Model information is contained in the set-up file along with 
other constants necessary in the data reduction process. There is a capability to do simple calcula- 
tions within the file which provides a mechanism for meeting unique computational test require- 
ments without having to implement special software. Reference 1 1 gives detailed information 
about the research data acquisition and reduction process. 

Off-line data reduction is done immediately upon receipt of the data reduction instruction 
sheet (Figure 9) from the test engineer or the research engineer. Tunnel customers and research 
engineers will be given a user area on one of the NTF file servers. The test data will be put into 
/tdata /ntf /testXX as soon as it is reduced. The customers/researchers may “ftp” the data files 
from the server to any workstation but they will not be given “write” privileges in this area. As 
previously mentioned, there is a two step process for data reduction in order to calculate the tun- 
nel flow angularity. If the value of the flow angle falls outside of the range 0.13 ±0.05 the 
researcher will be consulted. The data file name will have the form “ rxxx-yyyfp.f ’ where “xxx- 
yyy” is the run set,” f” indicates the presence of force data, “p” indicates the presence of pressure 
data and the “.f ’ appendix indicates that the flow angularity correction has been applied. The cus- 
tomer may have access to both corrected and uncorrected data if desired. Force and pressure data 
may also be separated into two files if this is necessary. Most research engineers prefer to do their 
own analysis but some post-run plotting is available at the customers request using the Data Anal- 
ysis System ( DAS) plotting tool (Reference 13) . This service is provided only if time permits 
and manpower is available. NTF test personnel will always monitor the health of the instrumen- 
tation and the quality of the data. Refer to section D for a discussion of data verification. 


B. Data Presentation 

The research data may be presented in many different ways. The Graphical User Interface 
(GUI) software (SAMMI ) allows for creating schematics and tabular displays. Several different 
types of plots and menus can be created for easy user interface to these displays. A real-time 
“snapshot” of selected raw and computed data will be configured pre-test and displayed on the 
workstation screen or printed to the lineprinter. Cyclic or “real-time” data is presented with the 
GUI but point-based data is usually displayed through the Research Graphics System (RGS). 

This system consists of third-party software (RT-WORKS) which provides for both continuous 
and point data to be plotted and displayed on as many as three Sun workstations. An additional 
feature of this system allows for previous test data to be plotted along with the current test infor- 
mation. This feature is very helpful in diagnosing problems and also in making “real-time” deci- 
sions about the test program. Once the RGS is configured for a test, the research engineer and/or 
customer will be allowed to change screens, change scales on individual plots and display pre- 
defined comparison data during operations. As with most on-line graphics systems, there are cer- 
tain restrictions. 

C. Data Transmittal 
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The information to be transmitted during the test will be determined pre-test and, if possible, 
will remain the same during the test. It is important to include all information which might need 
to be analyzed. The transmittal format is one that has been used at NTF for many years and is 
based on the standard interface file (SIFT) concept . The format is an ASCII 80 character record 
which can easily be read on any computer. The data items included in these transmittal files can 
be changed at any time during the test at the request of the customer through the proper channels. 

Each NTFSIFT file contains at least one set of NAME records. Following these records, there 
may be many sets of DATA records. There may or may not be other sets of NAME records on the 
file. 

NAME records are written with the following FORTRAN statements: 

WRITE (IOUT501) (ARRAY) I),I=J1,J2) 

501 FORMAT (2X,6(A8,5X)) 

where 


ARRAY = array of names 

J2 =J1 + 5 and J1 is incremented by 6 until J2 > number of names. 
All unused fields are filled with the name “BLANK " . 

DATA records are written with the following FORTRAN statements: 

WRITE (10UT, 502) (ARRAY(I),.I=J1, J2) 

502 FORMAT (2X, 6(E12.6, IX)) 

where 


ARRAY = array of data 

J2 = J1 + 5 and J1 is incremented by 6 until J2 > number of data items. 

All unused fields are set to -9990.0. 

The final data can be transmitted on 4mm DAT tape, optical disk or electronically transferred. 
The latter is the preferred mechanism. Flow visualization data may also be transmitted on CD- 
Rom. All research data generated in the NTF is FEDD restricted and company proprietary. 


D. Data Verification 

Because data quality is very important in the National Transonic Facility, data verification 
becomes the task of every individual involved in the test process. The data verification proce- 
dures take place not only when a model is being tested in the tunnel but when the model is placed 
in one of the model preparation areas. Data verification can logically be divided into two catago- 
ries: 
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(1) On-line Verification Procedures 

Online Snapshot 

A list of items will be predetermined for each test and these values will be saved in a 
snapshot file and/or printed on the lineprinter. These items are usually a combination of raw data, 
engineering units and computed data. Since a full set of corrections can be and usually are 
applied to the data, the values can be compared to the offline values post-run. The wind-off-zero 
point and tare run used to correct the balance data are also displayed. The snapshot is used to 
monitor the on-line calculations to insure instrumentation and data system health. 

Wind-off Zero Log 

A record of balance raw voltages, model angles calculated from the model acceleromenter 
and the arcsector accelerometer, balance temperatures, and tunnel temperature are recorded for 
each wind-off zero point. At the present time, this is done by the data system engineer but future 
plans are to have this stored on a computer file which can be retrieved at any time by test person- 
nel. This log monitors the health of both the balance and model orientation devices. 

Diagnostic Messages 

For purposes of maintaining the health of various tunnel sytems such as the ESP 8400 system, 
the NEFF DAU, and the MACH Ruska system, many diagnostic messages are printed either to 
the Data System Engineer (DSE) workstation or to the lineprinter to indicate a potential problem. 

Daily Instrumentation Calibrations/Checks 

Calibrations and checks are run each morning by the instrumentation technicians on all of the 
above mentioned systems. Calibration reports are checked and the test engineer is notified if any 
of the systems are out of tolerance. If hardware maintenance is needed, it can be done early 
before testing begins. A report of this activity is given at the morning shift meeting. 


OADTSP tool 

This tool provides a dynamic display of real-time channel data including the raw value, EU 
value, channel name, description, Neff range and filter setting. 

RAWDISP tool 

This tool can help provide verification of input data as close to the hardware as possible by dis- 
playing raw voltages (p V). 

On-line Research Data Plots 

Research plots are set-up pre-test by the research engineer and the test engineer. These plots 
can be a combination of real-time data and point-based data. The requested data is sent across the 
network from the RCS computer to a Sun workstation. There are three screens available for dis- 
plays and plots. The RGS allows for historical comparisons (test-to-test), repeatability compari- 
sons (run-to-run within a test), and logical progression checks. By viewing these “real-time” 
comparison plots, problems can often be found early enough to save nitrogen consumption. 
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RCS Debug Report 

If a problem is suspected, various debug reports can be turned on by the DSE. Many of the 
intermediate calculations are saved for these reports and problems can sometimes be resolved by 
displaying these intermediate values. The user should be warned that only those reports that are 
needed should be initiated because large files and/or large stacks of paper can be created through 
this process. 


(2) Off-line Verification Procedures 

Off-line verification of data involves running one or more of the following “check” programs 
available on the ModComp or the Sun workstations : 

RAWANAL 

This software is available on the ModComp and can be run by the DSE at the request of the 
test engineer, research engineer, or other test personnel. The raw data file is read and individual 
scans of data for a requested point or points can be displayed as well as the average raw data 
record. Either ESP or Neff data can be displayed. 

RAWED 

This software is available on the Sun workstation and can be used by the DSE to modify infor- 
mation in the raw data header records such as dataid, run and point. For example, records can be 
skipped or rearranged to eliminate a bad ESP calibration. 

DAUSCAN 

This software is available on the Sun workstation and can also be run by the DSE at the request 
of the test personnel. ESP or DAU channel or channels may be requested and the individual 
scans are displayed. Statistical information is computed for a particular channel including mean, 
standard deviation and other information that might be helpful in assessing data problems. This 
program is particularly useful in detecting noise and lag time problems in the data. 

TUNPC 

This software is available on the Sun workstation and can be used by any of the test personnel 
to check tunnel parameters. With minimal instruction from the data quality engineer, this tool is 
useful in checking the continued veracity of the thermodynamic and equation of state subroutines. 
This program has an option to correct the total pressure for height difference. 

TEMPCOMP 

This software is available on the Sun workstation and can be used by any of the test personnel 
to check temperature compensation. First, the temperature compensation file for this test must be 
downloaded from RCS to the workstation. The program calculates temperature compensation for 
the six balance components for any given data point. The wind-off zero data for that point is also 
necessary. 
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This software is available on the Sun workstation and can be used by any of the test personnel 
to check 1st and 2nd order balance interactions. First, the latest balance calibration interaction file 
must be downloaded from RCS to the workstation. Then the balance forces and moments cor- 
rected for interactions can be calculated. 

EULER 1 

This software is available on the Sun workstation for all test personnel. Given input angles in 
the gravity axis system and rotation information, it provides a check on angle rotations from the 
gravity axis system to the balance axis system and then the model axis system. 

EULER2 

This software is available on the Sun workstation for all test personnel. Given forces, moments 
and angles in the gravity axis system and the rotation information, it provides a check on forces 
and moments in the body axis, the stability axis and the wind axis systems. 

The NTF operations personnel will be happy to assist the tunnel users in data verification using 
any of the above mentioned procedures. The ultimate product of the National Transonic Facility 
is high quality research data. The methods for acquiring, reducing and presenting research data in 
the NTF are highly flexible. Every effort will be made by the NTF staff to honor the requests of 
our customers when time and manpower are available. All data in the National Transonic Facility 
are considered FEDD restricted and company proprietary and certain restrictions are necessary 
to protect the data. 

E. Future Plans 

Future plans for the data system include: 

1. On-line flow angularity calculation and corrections 

2. On-line wall corrections 

3. On-line research data used by model attitude control algorithm to achieve desired 
angle of attack (Alpha)and angle of sideslip (Beta) 

4. On-line calculation, display and recording of standard deviation for all Neff channels 

5. Off-line relational database storage and retrieval of test data for specific programs. 

vra. CONCLUSION 

In this document, the users of the NTF are provided with information about the data system, 
the various instrumentation devices, the data reduction algorithms, the data management system 
and the data transmittal procedures. This report only briefly describes the design and major char- 
acteristics of the tunnel, and it is anticipated that it will be used in conjunction with other docu- 
ments that fully describe these aspects of the National Transonic Facility. A new version of the 
National Transonic Facility User’s Guide can be found at the following address on the World 
Wide Web: 


http://ad-www.larc.nasa.gov/facility/ntf/ntftb.html 
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APPENDIX A 


BALANCE DATA REDUCTION EQUATIONS 


Three types of data are required to properly compute the loads measured by a strain-gage bal- 
ance: wind-off reference data, weight tare data, and wind-on data. Wind-on is the term used to 
refer to data acquired to measure an unknown load due to forces other than the model weight. 
The processing used to compute the unknown load for wind-on data requires the prior processing 
of the wind-off and the weight tare data. A detailed mathematical explanation of this procedure 
will follow. 

The processing of a wind-off data point involves the calculation of a set of weight tare factors 
which are used to compute the initial loads on each balance. These factors are calculated by an 
iterative procedure. A mathematical explanation of this procedure is included in the following 
paragraphs. 


The model weight is a vector, W . In the gravity axis system, it is apparent that there are no hor- 
izontal components of weight and that the vertical component is directed downward. The weight 

vector in the gravity axis system, , is then 

W s = (W gx ,W gy ,W gz )= (0,0, -WO 


where W is the magnitude of W s at attitude [ R gb ] . This attitude matrix is defined as follows: 

IV- [**♦,] IW !• w 

where <J> g , a g , y^ are rotations about the roll, pitch and yaw axes respectively and the subscript 
denotes the gravity axis system where: 




(cosa^cosyp 

(cos<}> siny -sin<|) sina cosy ) 
8 8 8 8 8 

(sin<]) sin\|/ + cos<J) sina cos\|/ ) 

o <5 o o o 


(-cosa^sinyp 

( cos <|> cosy + sin<(> sina sin y ) 
8 8 8 8 8 

( sin cos y g - cos ^ sin a g sin y^ ) 


(-sinap 
(-sin^cosap 
(cos<|> cosa ) 

O O 
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Because the x and y components of W g are zero, only the third column of [/? 
sidered. 


Force Weight Tare Computation 

Consider the third column of [R gb ] , denoted [R3] , which is 


[/?3] = 


-sina 

8 

-sin<t> cosa 
T 8 8 

cos (f> cos a 

- 6 oj 


For the wind-off zero point [/?3] becomes [/?3 0 ] : 


[*3 0 ] = 


-sina 0 

-sin<() 0 cosa 0 

cos<|> 0 cosa 0 


The delta weight tares are defined as the change in balance components due 
W, relative to initial loads. Now let 


[H= -([/?3] - [K3 0 ]) 

Then the delta force weight tares can be expressed as: 


AFT A 


VI 


W ( sina^ - sina 0 ) 

SFTA 

= W 

V2 

= 

W (sin^coso^ - sin<(> 0 cosa 0 ) 

NFTA_ 


V3_ 


-W (cos^cosa^- cos<|> 0 cosa 0 ) 


Now, recognizing that the force beams of a balance will each sense a differc 
balance weight as well as the model weight, this equation may be rewritten as: 


WAF= 


AFT A 
VI 


WSF= 


SFTA 

V2 


] needs to be con- 


model weight 


portion of the 
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WNF= NFTA 

V3 


provided the denominators are not zero. 

In actual practice, several weight tare recordings are made at different balance attitudes which 
results in an overly defined system of equations. The solution for the force weight tares will be 
written in the following form for consistency with the moment weight tares. 

Assuming n weight tare recordings, the solution is written: 


WAF= 


X (AFTA'XVl') 


i = i 

DTAF 


n 


WSF= 


X (SFTAj x V2 { ) 

i = 1 


DTSF 


where the denominators are: 


WNF= 


X {NFTA i x V3.) 


is 1 

DTNF 


n 

DTAF= 2 (Vl ( .x VI.) 

is 1 


n 

DTSF= ^ ( V2 i x V2 t ) 

i = 1 


n 

DTNF= J ( V3 f x V3 ( .) 

i = i 


Moment Weight Tares 

If transfer distances X, 7, and 2 are measured in the balance axis system from the balance 
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moment center to the model center of gravity, positive in the directions of positive thrust, side 
force and normal force respectively, then the delta moment weight tares are obtained by transfer- 
ring moments as follows: 


RMTA 


SFTAxZ-NFTAxY 

PMTA 

— 

AFTA xZ + NFTA xX 

YMTA_ 


SFTA xX + AFTA xY 


Substituting for AFTA, SFTA and NFTA gives: 


RMTA 


(WSFx V2 x Z) - (WNFxV3xY) 

PMTA 

~ 

( WAF x V\ x Z) + ( WNF x V3 x X) 

YMTA_ 


JWSFxV2xX) + (WAFx VI xY)_ 


Each of the equations can be solved for the weight tare factors by first premultiplying both sides 
of the equation by the transpose of the row vector on the right containing the V terms. The solu- 
tion then involves inverting the matrix obtained from the transpose of the row vector times the 
row vector itself. 


Recalling that in actual practice, several weight tare recordings are made (both pitch and roll 
angles), the solution obtained are then: 


X ( V3 i x V3 ,) x X ( RMTA i x V2 ,) “ X ( V2 i xV3 i) x X ( R MTAjXV 3 y ) 


WZRM= 


i = 1 


/ = I 


i = 1 

DTRM 


/= 1 


WYRM= 


X ( V2 i xV3 i) x X ( R MTA f xV 2.) - £ (V2 i xV2 j ) x £ (RMTA j x V3 ( ) 


/ = 1 


7 = I 


7 = 1 


7 = I 


Z)77?M 


WZPM= 


X ( V3 / x V3 ,) x X (PMTA.xVI.) - £ (V1 ( .X V3 ( .) x £ (PMTA.xV. 3 ( ) 


7 = 1 


1= I 


7 = 1 


DTPM 


WXPM= 


n 


n 


X (W/Xl'l,) x X ( PMTA i xV3 i ) - £ (Fl ( .x 

LzJ LzJ )' = 1 

DTPM 


V3 .) x £ (PMTA t x VI.) 

i=l 
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X ( VI . x VI,.) x £ ( FA/771, x V2,) - £ ( V2 i x VI ,) x £ ( YMTA i x VI ,.) 
WXYM= ^=2 


/ = 1 


/ = 1 
DTYM 


i = ] 


X ( V2 , x V2 ,) x X ( x V] ,) - X ( V2 - x V1 ,) x X ( YMTA > x V2 /) 

WYYM= 


/ = I 


i = 1 

DTYM 


i = 1 


where the denominators are 


DTRM= J^(V2 x V2,) x £ (V3,x V3 f ) - X (V2 f x V3 f ) x £ (V2, x V3.) 


f= 1 


I a 1 


0™= X ( V1 , XV1 /) x X ( V3 i xV3 i) - X ( vl i xV,3 ( ) x X (Vl.-x V3-) 


I a 1 


i = 1 


I a 1 


f = 1 


DTYM= X (Vl f x VI,) x X (V2,-x V2 f ) - X ( Vl , x V2 ( ) x £ (VI, x V2,) 


i = i 


i = I 


i = I 


i = 1 


If the absolute value of the demoninator for either the force or moment equations is less than 
some predetermined tolerance level, then this is an indication that there is insufficient data in this 
set of tare points to accurately calculate this particular parameter. The parameter is then set equal 
to a corresponding parameter for which there is sufficient information to be calculated (Reference 
12 ). 


The nine weight tare factors are computed based on corrected balance-axis components, which 
in turn are based upon assumed initial loads equal to the wind-off zero loads. The newly com- 
puted weight tare factors are then used to recompute the initial loads, which are then compared to 
the original loads. If the new initial loads are sufficiently close to the old initial loads, they are 
assumed to have converged, and the tare computations are accepted. If the new initial loads are 
significantly different, the old initial loads are replaced with the new initial loads, which are then 
used to recompute the corrected balance components to start another iteration. This process con- 
tinues until the initial loads converge or for a maximum of five iterations. The accuracy required 
for each component is obtained from the balance calibration. The nine weight tare factors, the six 
initial loads, and the wind-off zero angles are all saved for future wind-on computations. 


Corrections for Balance Interactions 

The uncorrected balance components are AF, SF, NF, RM, PM, and YM. These components 
are referred to as “delta” components because they are relative to a wind-off zero recording of ini- 
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tial loads. The components are treated as a 6 X 1 column vector demoted by FU (forces uncor- 
rected), where: 

AF 

SF 

FU= NF 

RM 

PM 

YM 

Let the components corrected for interactions be denoted [ F ] 


[F] = 


AFC 

SFC 

NFC 

RMC 

PMC 

YMC 


For first order interactions, the following matrix relationship between correct and indicated 
components has been established during the balance calibration: 


[FU]= [Cl] [F] 


where [Cl ] is a 6 X 6 matrix which is the normalized first order interaction coefficient matrix 
with main diagonal elements of unity. 

Provided [Cl] is nonsingular, the correct delta components are found from: 


[F]= [Cl/] [FU] 


where [Cl/] is the inverse of [Cl]. 

Second Order Interactions 


For second order interactions, the following matrix relationship between correct and indicated 
components has been established during the balance calibration: 

[FU]= [Cl] [F] + [C2] [F2] 
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were [ C 2] is a 6 X 2 1 matrix which is the normalized second order interaction coefficient matrix 
and [F 2] is a 21 X 1 matrix of product combinations of [F] as follows: 


[F 2] = 


AFCxAFC 
AFCxSFC 
AFC x NFC 
AFCxRMC 
AFC x PMC 
AFCxYMC 
SFC x SFC 
SFC x NFC 
SFCxRMC 
SFC x PMC 
SFC x YMC 
NFC x NFC 
NFC x RMC 
NFC x PMC 
NFC x YMC 
RMC x RMC 
RMC x PMC 
RMC x YMC 
PMC x PMC 
PMC x YMC 
YMC x YMC 


Provided [Cl] is nonsingular, The correct delta components are found from: 

[F]= [C\l] [FU] - [C1/C2] [F2] 

where [C1/C2] is the product of the inverse of [Cl] and [C2] . This equation must be solved 
iteratively because [F] is expressed in terms of [F2] which is itself a function of [F] . 

Translation of Initial Loads 

Typical balance calibration coefficients as determined in the laboratory have a zero-load refer- 
ence (i.e., no external weight or load applied to the balance) and, therefore, are based on total 
applied loads. In the wind tunnel, however, load measurements are made relative to a wind-off 
zero recording. These loads are delta loads since the weight of the model has already been sensed 
by the balance when the wind-off zero was taken. Since the 2nd-order interaction coefficients are 
nonlinear it is necessary to add in the initial loads before correcting for 2nd-order interactions. 
The following is a mathematical description of the procedure. 


37 



Let [/ r i/ 0 ] denote the indicated initial loads, 
[F 0 ] denote the correct initial loads, 

[FUT\ denote the indicated total loads, and 
[FT] denote the correct total loads 

where 


[FT] = 


AFTO 
SFTO 
NFTO 
RMT 0 
PMT 0 
YMTO 


By definition 


[FUT}= [FU] + [F 0 ] 
and 


[FT}= [F] + [F 0 ] 

The second order interaction relationship between correct and indicated components for the tun- 
nel is then: 


[FT]= [Cl/] [FUT] - [CIIC2] [F2T] 

where [F2T] is just [F2] based on [FT] . This equation must be solved iteratively. 
Note that for initial loads, this becomes 

[F 0 ]= [Cl/] [Ft/ 0 ] - [C1/C2] [F2 0 ] 

where [F2 0 ] is just [F2] based on [F 0 ] . 

Now, for the i ,h iteration, let 


[£;]= [C1/C2] [F2T] 


and specifically, for initial loads, let 
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[e 0 ]= [C1/C2] [F 2 0 ] . 


The iteration technique is then given by: 


Iteration 

1 FT] Approximation 

Error 


0 

[CM] [ FU ] + [CM] [FU 0 ] 

[e 0 ] 


1 

[CM] [FU] + [CM] [FU 0 ] - [e 0 ] 

[£,] - 

[e 0 ] 

2 

[CM] [FU] + [CM] [FU 0 ] - [e,] 

[e 2 ] - 

[e,] 


This iteration is continued until [e ; ] is less than the specified accuracy for all compo- 

nents. The accuracy criteria are obtained for the balance calibration. 

The correct delta components are then obtained from; 

[F]= [FT] - [F 0 ] 


Computation of Sting Deflections 

Sting deflection (or bending) occurs when loads are applied through the model. Devices used 
to measure angles which define the balance and model attitude are sometimes located such that 
they do not record deflections of the sting. Consequently, sting deflection angles must be com- 
puted as a function of correct balance loads. The deflection angles are assumed to be small 
enough so that the sting responds elastically, allowing the angle to be described as a spring con- 
stant times a load. The loads used in the bending calculations can be either total loads or delta 
loads. In the following equations the spring (deflection) constants are SFDF, YMDF, NFDF, 
PMDF, and RMDF. 

For delta loads, the current deflections and initial deflections are: 

YA WS= (SFC x SFDF ) + ( YMC x YMDF) 

ALPS= (NFC x NFDF) + {PMC x PMDF) 

PHIS= ( RMCxRMDF ) 

YA Wl- 0.0 
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ALPI= 0.0 


PHII= 0.0 

For total loads, the current and initial deflections are: 

YA WS= (SFTO x SFDF ) + ( FA/7’0 x YMDF) 
ALPS= (NFTOxNFDF) + (PMTOxPMDF) 
PHIS= (RMTOxRMDF) 

YAWI= - YAWS 
ALP!= -ALPS 
PHII= -PHIS 


Computation of Balance Attitude 

The attitude of the balance with respect to gravity is determined on the basis of a specified 
input rotation scheme which consists of an ordered set of orthogonal Eulerian transformations 
(Appendix B). Each transformation rotates the components of a vector through a specified angle 
about a specified axis. The final result is the transformation of the components of a vector from 
the gravity axis system to the balance axis system. The axis system used are right-hand Cartesian 
systems (Figure 10). 

The transformation matrix , as defined on the first page of this appendix, describes the 

attitude of the balance with respect to gravity and is derived from successive applications of the 
appropriate individual transformations. This can be summarized as a single yaw rotation YAWG, 
followed be a single pitch rotation ALPG, followed by a single roll rotation PHIG. That is: 

[/?**]= [R X (PHIG)] [ R y (ALPG )] [/? z (L4WG)] 

Solving for YAWG, ALPG, and PHIG gives: 

YAWG= atan [ (cosa sin vp )/(cosa cos\|/ )] 

6 6 & & 


ALPG = asin [sina 1 
L s J 

PHIG= atan [ (sincf) cosa )/ (cos <)) cosa )] 
s s s s 
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Restrictions and special conditions for these calculations are listed in Reference 12. In the NTF 
data reduction software there can be up to eight gravity to balance rotations with a special code to 
define the type of rotation. This allows for rotations due to a bent sting (knuckle). 

The gravity to balance matrix [R gb ] is then further rotated through the sting-bending angles 
PHIS, ALPS, and YAWS by the method described in APPENDIX B. 


Correction for Weight Tares 

For the balance at any attitude [K , let [FT ARE] equal a vector containing weight tare cor- 
rections as defined previously. 


[FT ARE] = 


AFT A 

SFTA 

NFTA 

RMTA 

PMTA 

YMTA 


Now let [FBAL] denote the aerodynamic loads in the balance axis system where: 


[FBAL] = 


AFBA 

SFBA 

NFBA 

RMBA 

PMBA 

YMBA 


The aerodynamic loads in the balance axis system are computed by subtracting the weight tares 
from the correct delta balance loads. 


[FBAL] = [F] - [FT ARE] 


Computation of Model Attitude 


Now that the gravity axis to balance axis transformation has been fully developed. 

Similar transformation matrices from balance axis to model axis [ Rj )m ] and from wind axis to 
gravity axis [R ] may be similarly defined on the basis of specified input rotation schemes. 
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The wind to gravity matrix [R wg ] describes the tunnel flow angularity and is defined as an 
upflow or pitch rotation and a crossflow or yaw rotation. 


The attitude of the model with respect to the wind can then be described by the transformation 
matrix given by: 


i R w J= [**J tv tv 

The model attitude [/? vvm ] can also be summarized as a single roll rotation PHIW, followed by 
a single yaw rotation YAWW, followed by a single pitch rotation ALPW. That is: 

tV= [R y (ALPW)] [/? z (L4WW)] [RJPHIW)] 

Substituting the elementary transforms and carrying out the indicated multiplications yields: 


R 


WM 


(cosa„cosv| /J 
(sinyj 
(sina^cosij/J 


(-cos a M sin cos \\i w - sin a w sin<|> J 
(cos^cos^J 

(-sina M ,sinv|^cos<t> w + cosa w sin<|>J 


( cos a w sin sin <|> H , - sin a w cos cp J 
(-cos\]/ H ,sin<t>J 

( sin a w sin sin <j) ^ + cosa^cos^,) 


Solving for YAWW, ALPW, and PHIW gives: 

PH1W= atan [ (cosy^sin^) / (cosxj/^cos^) ] 

YAWW= asin [sin\j/ M ,] 

ALPW= atan [ ( sin ct w cos \\i w ) / ( cos a w cos \\i w ) ] 

The restrictions and special conditions associated with these angle computations are listed in Ref- 
erence 12. 

The angle of attack of the model is ALPW and the angle of sideslip is computed from the defini- 
tion (Figure 11): 


BETA= -YAWW 


42 



Computation of Balance Components in the Model Axis System 


The correct balance axis components can now be rotated into the model axis system by using 
rotation matrix [R bm ] . The subscripts are used for clarity in describing the final rotations and 
corrections to the balance components. 


AFMA ] 


AFBA 

SFMA ] 

= 

SFBA 

NFMA X 


NFBA_ 


-RMMA , 


-RMBA 

PMMA X 

fl 

>3 
3 ^ 

PMMA 

-YMMA , 


-YMBA 


The minus signs in the moment equations result from axis system differences in the conventional 
wind tunnel balance force and moment sign conventions. 

The transfer distances XBAR, YBAR, and ZBAR describe the transfer from the balance moment 
center to the model moment reference center (Figure 12). The transfer equations are given by: 

AFMA 2 = AFMA { 

sfma 2 = sfma 2 
nfma 2 = nfma x 

RMMA 2 = RMMA , + ( NFMA X x YBAR) - (SFMA l x ZBAR) 

PMMA 2 = PMMA { - ( NFMA } xXBAR ) - (AFMA^xZBAR) 

YMMA 2 = YMMA { - ( SFMA { xXBAR ) - ( AFMA X x YBAR) 

Base and Chamber Pressures 

Letting PB i denote the i h base pressure, the base pressure coefficients are computed as: 
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where 


and 


q m = free-stream dynamic pressure 
= free-stream static pressure. 


The forces and moments to be applied as a correction are computed as: 


AFB 

SFB 

NFB 

RMB 

PMB 

YMB 


NBAS 

£ ( PB-PJ 

i = I 


-AREA.g. 

AF , 

area sf 

area nf 

(AREA x ARM) RM 
(AREA x ARM) PM 
(AREA xARM) ym 


where the summation i is over all base pressures for which the correction flag is set equal to one. 
The column vector on the right side and the correction flag are obtained from the input specifica- 
tions. 


If the correction flag is set to zero, the forces and moments (XAFB, XSFB, XNFB, XRMB, 
XPMB, XYMB) are computed as above but not applied. 

The axial force terms are also expressed in terms of axial and drag coefficients as: 

CAB= AFB/(q oo xS) 


CDB= CAB x cosa 
XCAB= XAFB/ ( g<x> xS) 


XCDB= XCAB x cosa 

W 


Letting PC i denote the i ,h 


chamber pressure, the individual pressure coefficients are computed 



CPC,= 


( 'PC-PJ 


The forces and moments to be applied as a correction are computed as: 


AREA af 
area sf 

AREA nf 

(AREA x ARM) RM 

(AREA xARM) pM 

(AREAxARM) ym 

where the summation of i over all chamber pressures for which the correction flag is equal to one. 
The column vector on the right side and the correction flag are obtained from the input specifica- 
tions. 

When the correction flag is equal to zero, the forces and moments (XAFC, XSFC, XNFC, 
XRMC, XPMC, and XYMC) are computed as above but not applied. 


AFCH 

SFCH 

NFCH 

RMCHl 

PMCH 

YMCH 


NCHM 

- I <pc,-p _) 

i = 1 


Sign Convention for Area (A) and Arm (1) in Base and Chamber Correction 


Force/ 

Orifice 

Area (A) 

Arm (1) 

A*I 

Moment 

Location 

Sign 

sign 

Axial 

Rearward facing 

+ 



Axial 

Forward facing 

- 



Normal 

Top 

+ 



Normal 

Bottom 

- 



Pitch 

Top-rear 

+ 

- 

- 

Pitch 

Bottom-rear 

- 

- 

+ 

Pitch 

Top-forward 

+ 

+ 

+ 

Pitch 

Bottom-forward 

- 

+ 

- 


The axial force terms are also expressed in terms of axial and drag coefficients as: 

CAC= AFCH/ (q„xS) 
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CDC= CACx cosa^ 
XCAC= = XAFC/(q„xS) 


XCDC= XCACx cos a 

W 


Computation of Model, Stability and Wind Axis Components 

The forces and moments are corrected for base and chamber pressures to give the final cor- 
rected model axis components: 


r 


ARM A, 


_ _ 



AFMA 


l 


AFB 


AFCH 

SFMA 


sfma 2 


SFB 


SFCH 

NFMA 


nfma 2 


NFB 


NFCH 

RMMA 


rmma 2 


RMB 


RMCH 

PMMA 


pmma 2 


PMB 


PMCH 

YMMA_ 


ymma 2 


YMB 


YMCH 


The stability axis components are obtained by rotating the model axis components through 
minus the angle of attack 


DRAG 

SFSA 

LIFT 




AFMA 

SFMA 

NFMA 



L 

J 

-RMSA 


-RMMA 

PMSA 

= [*,(-«*)] 

PMMA 

-YMSA_ 


-YMMA 


The wind axis components are obtained by rotating the stability axis components through the 
angle of sideslip: 


AFWA 

SFWA 

= [* Z (P)J 

DRAG 

SFSA 

NFWA_ 


LIFT 


-RMWA 

PMWA 

= [A (P)] 

-RMSA 

PMSA 

-YMWA_ 


-YMSA 
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Computation of Coefficients 


Let [FMA] , [FSA] , and [FWA] denote the model, stability and wind axis components and 
define the 6X6 maindiagonal matrix [ C] as: 


where 


[C] = 


1 


<7oo X5 


0 0 


< 7~ x5 


0 0 


0 0 
0 

0 0 


4oe X5 


0 


0 

0 

0 

1 


q xSxB 
0 0 0 


0 

0 

0 

0 

1 


0 


0 


q xSxCBAR 

1 oo 

0 


0 

0 

0 

0 

0 

1 


q^xSxB 


S 

CBAR 

B 


model reference area (ft ) 
model reference chord (in.) 
model reference span (in.) 
dynamic pressure (psf) 


and all of the elements off the main diagonal are zero as shown. 

Let [CM/4], [C5/4] , and [CWA] denote the model, stability and wind axis coefficients: 


[ CM A ] = 


CA 

CY 

CN 

CRM 

CM 

CYM 
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[CS/t] = 


CD 

CYS 

CL 

CRMS 

CMS 

CYMS 


[CWA] = 


CDW 

CYW 

CLW 

CRMW 

CMW 

CYMW 


The model, stability and wind axis coefficients are then computed as 

[CMA]= [C] [ FMA ] 


[CSA]= [C] [FSA] 
[CWA]= [ C][FWA ] 


The lift to drag ratio and lift squared are also computed. 


L/D= CL/CD 
CLSQ= CLxCL 
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APPENDIX B 


EULER ROTATION SCHEME 

The NTF data reduction software uses a rotation scheme to determine balance attitude, model 
attitude and also to adjust for flow angularity in the tunnel. This scheme (Figure 10) consists of 
an ordered set of orthogonal Eulerian transformations. Each transformation rotates the compo- 
nents of a vector through a specified angle about a specified axis. 

The following define the three types of rotations that are used: 

(1) Roll rotation through some angle about the X-axis - right wing down is positive. 

1 0 0 
Rx (ty) = o cos (<[>) -sin (<()) 

0 sin (<j>) cos (<()) 

(2) Pitch rotation through some angle 0 about the Y-axis - nose up is positive. 

cos (0) 0 -sin (0) 

Ry(9)= o i o 

sin (0) 0 cos(0) 

(3) Yaw rotation through some angle \\i about the Z - axis - right wing aft is positive. 

cos (\|/) -sin (y) 0 
Rz(y)= sin (\y) cos (v|/) 0 
0 0 1. 

The order of the rotations is important since the rotations are noncommutative. An ordered set 
of codes describing the axes of rotation is included in the data reduction set-up file. 
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APPENDIX C. 


DESCRIPTION OF NTF BALANCES 

Force balances for the National Transonic Facility must be specially designed to withstand the 
extreme temperature and high loads imposed by a cryogenic high-pressure tunnel. If a model or 
test requires a balance of a different size or load range then those listed in the table below, a min- 
imum of one year is needed to design, fabricate, gage and calibrate the balance for the NTF. 
However, if testing is to be done only at ambient temperature, the data reduction equations pro- 
vide for processing output from any type of balance, either force or moment. 


At this time all NTF balances are constructed of an 18-percent nickel maraging vacuum-remelt 
steel with an ultimate strength of 200 KSI. The model balance interface is a diameter fit with a 
dowel for fixing the longitudinal location of the moment center. The sting balance fit is a taper 
with set screws at 45° to seat and unseat the taper. In addition to the force and moment gages, 
each NTF cryogenic balance comes with three type T thermocouples positioned at the front, axial 
and rear of the balance. The output from these thermocouples is used to temperature compensa- 
tion the balance during testing. ( The NTF 1 1 3C has three PRT’s instead of the type T thermocou- 
ples. This balance also provides for wiring for the AOA package and ESP modules through the 
middle of the balance in order to keep these wires from carrying loads when they are passed over 
the outside of the balance. ) 

Most NTF balances have integral wiring 20 feet long. There are eighteen #32 gauge wires and 
three type T thermocouple wires encased in a sleeving. This wire bundle is approximately 0.25 
inches in diameter. The table on the next page lists all of the NTF balances with size and maxi- 
mum loads. 
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Size and Load Capacity Information on NTF Balances 


BALANCE 

SIZE 

diameter x length 

normal 

(lbs) 



m 

yaw 

(in-lbs) 

side 

(lbs) 

NTF101A 

2.375 X 15.565 

6500 

700 

13000 

9000 

6500 

4000 

NTF101B 

2.375 X 15.565 

6500 

700 

13000 

9000 

6500 

4000 

NTF 102 

2.000 X 14.281 

3000 

600 

6000 



300 

NTF 103 

2.000 X 14.281 

1500 



300 

300 

150 

NTF 104 A 

2.000 X 14.281 

3400 

300 

10000 

5000 

5000 

1000 

NTF104B 

2.000 X 14.281 

3400 

300 

10000 

5000 

5000 

1000 

NTF105 

2.000 X 14.281 

2000 

175 

6000 

3000 

3000 

700 

NTF106 

1.750 X 12.176 

2500 

350 


2500 

4000 

1000 

NTF 107 

0.750 X 7.250 

160 

50 

250 

100 

125 

40 

NTF 108 

1.500 X 12.013 

1600 

125 

3000 



500 

NTF 109 

Nose Balance 

1200 

175 



1500 

350 

NTF111 

1.750 X 12.176 

1000 

BB 


1000 

2000 

500 

NTF112 

1.250 X 9.400 

600 

180 

1700 

1000 

400 

280 

NTF113A 

2.375 X 15.565 

6500 

400 

13000 

9000 

6500 

4000 

NTF113B 

2.375 X 15.565 

6500 

400 

13000 

9000 

6500 

4000 

NTF113C 

2.375 X 15.565 

6500 

400 

13000 

9000 

6500 

4000 

NTF 11 3D 

2.375 X 15.565 

6500 

400 

13000 

9000 

6500 

4000 

NTF114S* 

16.00X25.750 

6100 

1300 

70000 

353000 

75400 

l 

i 

i 

i 

NTF115 

2.000 X 14.281 

3200 

200 

B1 


2500 

1000 


* The NTF1 14S balance is a side-wall balance used for semi-span testing. 








































































































































APPENDIX D. 


ESP 8400 SYSTEM DIAGNOSTICS 

The following diagnostic messages are part of the standard RCS ESP calibration software. 
The requirements apply to full module calibration unless otherwise indicated. (FS = full scale 
pressure of the module and C are the calibration coefficients.) 

test id 

|C 0 | < 0.25 X full scale pressure of module 

MESSAGE m 

C 0 out of tolerance ( C 0 > 0.25 X FS ) 

TEST (21 (zero calibration) 

|C 0 | < 0.004 X FS 

MESSAGE (2) 

AC 0 out of tolerance ( C 0 > 0.004 X FS ) 

TEST (3) 


[C,> (FS/5) x0.7] .AND. [ (C, < FS/5) x 1.25] 

MESSAGE (31 (if out of these bounds) 

C, (sensitivity) out of tolerance 

TEST (4) 

The combined non-linear coefficient test should be based on the maximum calibration pressure 
voltage for the given port (i.e. V4). 


C 2 x (V4) 2 + C 3 x ( V4) 3 x C 4 x (V4)‘ 


<0.007 X FS 


MESSAGE (4) 

Combined coefficient calculation out of tolerance 
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Figure 1-A. National transonic facility 
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Figure 1-B. Plan of tunnel circuit. 
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Figure 2- A. Tunnel Operating Envelopes 
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Figure 2-B. Tbnnel Operating Envelopes (continued) 
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Figure 2-C. T\innel Operating Envelopes (continued) 
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Figure 2-D. Tunnel Operating Envelopes (continued) 
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Figure 2-E. Tunnel Operating Envelopes (continued) 













NTF ENVELOPE~MACH=0.600 



VISd ‘3UnSS3UdlV±0± 


62 


Figure 2-F. Ttonnel Operating Envelopes (continued) 
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Figure 2-G. l\innel Operating Envelopes (continued) 



















NTF ENVELOPE— MACH=0.800 
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Figure 2-H. Ttinnel Operating Envelopes (continued) 
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Figure 2-1. Tunnel Operating Envelopes (continued) 
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Figure 2-J. Ttinnel Operating Envelopes (continued) 
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Figure 2-K. l\innel Operating Envelopes (continued) 
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Figure 3. Model access system. 
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Figure 4-A. Research Computer System Hardware Schematic 





Display Workstations 
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Figure 4-B. Process Computer System Hardware Schematic 
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Figure 4-C. Model Preparation Area Computer System Hardware Schematic 




















PCS Block Diagram 
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Figure 5-A. Process Computer System Functionality Diagram 








74 


Figure 5-B. Research Computer System Functionality Diagram 















MODEL ACCELEROMETER - CASE 



1 . BOLONCE IN MODEL OX IS SYSTEM. 0 BM '= 0 

2. BOLONCE UPRIGHT, (f) BM = 0 


Figure 6-A. Model Orientation Definitions using Model Accelerometer 
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MODEL ACCELEROMETER - COSE TL 



1 * BALANCE in model ox is system, 0 = o 

BM 

2. BALANCE INVERTED, (f> = ] 80 


Figure 6-B. Model Orientation Definitions using Model Accelerometer 

(continued) 
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MODEL ACCELEROMETER - COSE JR 



1. MODEL ACCELEROMETER MISALIGNED UITH RESPECT 

TO MODEL AXIS, Q QMS (GRAVITY AXIS # MODEL AXIS) 

2. BALANCE UPRIGHT, (f) BM = 0 


3. BALANCE MISALIGNED UITH RESPECT TO MODEL REFERENCE 
AXIS, Q RM (BALANCE AXIS f MODEL AXIS) Kt= ^ KLNLt 


Figure 6-C. Model Orientation Definitions using Model Accelerometer 

(continued) 
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MODEL ACCELEROMETER - COSE ZESL 



1. MODEL ACCELEROMETER MISALIGNED WITH RESPECT 

TO MODEL AXIS, Q 0MS (GRAVITY AXIS * MODEL AXIS ) 

2. BALANCE INVERTED, (f) BM = 180 


5 . 


BALANCE MISALIGNED WITH RESPECT TO MODEL 
AXIS. Q bm (BALANCE AXIS # MODEL AXIS) 


REFERENCE 


Figure 6-D. Model Orientation Definitions using Model Accelerometer 

(continued) 
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GRCSECTOR CfiSE JE 


1. STING OFFSET, 9 k 

2. BALANCE UPRIGHT, (j) BM = 0 



9 M0D ~ &BRCSEC + ®k ®SB 
®MOD = ®ACC~ ®SM + ®k~ ®SB 


Figure 6-E. Model Orientation Definitions using Arcsector Accelerometer 
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ARCSECTOR CfiSE YL 


1. STING OFFSET, 0 k 

2. BALANCE UPRIGHT, (j) = 0 

_ ' BM 

3. BALANCE OFFSET FROM MODEL REFERENCE, Q RM 

4 . STING BENDING, e sB 



Figure 6-F. Model Orientation Definitions using Arcsector Accelerometer 

(continued) 
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ARCSECTOR CASE JEZT 


1. 

2 . 

3 . 

4 . 


STING OFFSET. 0 k 

BALANCE INVERTED. ^ BM = 1 80 ' 

BALANCE OFFSET FROM MODEL REFERENCE. 

STING BENDING. Q sb 


0 


BM 




Figure 6-G. Model Orientation Definitions using Arcsector Accelerometer 

(continued) 
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COMBINED ANGLE SYSTEMS - CASE 1 THL 


1. MODEL ACCELEROMETER MISALIGNED WITH 
RESPECT TO MODEL AXIS, Q 0MS 

2. BALANCE INVERTED, (j) =180' 

3. BALANCE MISALIGNED WITH RESPECT TO MODEL 
REFERENCE AXIS, 0 OA , 

BM 



Figure 6-H. Model Orientation Definitions using Arcsector Accelerometer 

(continued) 


82 











83 



National Transonic Facility 
Test Section Model 
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Figure 8. ESP System Schematic 
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Figure 9. Sample Data Reduction Instruction Sheet 
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Figure 12. Transfer distance definitions for balance to model moment center translations. 
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